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H , Abstract 



We perform a 1/A^-expansion of the partition function of the massive Gross- Neveu 
model in 1+1 dimensions. The procedure allows for the inclusion of the contribution 
of scalar and pseudoscalar composites (of order 1/A^) to the equation of state. The 
naive expectation that the bosonic fluctuations correct significantly the mean field 
approximation at low temperatures is confirmed by our calculations. Actually the 
relevant degrees of freedom of hadronic matter at low temperatures are found to be 
pion-like excitations, rather than the fundamental constituents. 
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1. Introduction 



The equilibrium thermodynamics of strong-interacting matter has been studied, up 
to now, by quite different approaches. Apart from lattice simulations, most of the other 
studies rely on phenomenological models, which try to capture some relevant physical 
features of the theory |^, ^ . The absence of a common starting point in the literature 
is due to the non perturbative character of low energy QCD, and to the vicinity of the 
expected phase transition, at temperatures of the order of Aqcd- Thus, in employing 
analytical methods to study the transition, one is faced by this major difficulty of QCD, 
that of passing from the hadronic regime to that of quarks and gluons, based on the same 
fundamental lagrangian. 

For very low temperatures, important results have been obtained on the assumption 
that the thermodynamics of hadronic matter in this regime is that of a gas of weakly 
interacting pions [§], since pions are the lowest mass excitations and are thus expected to 
dominate the partition function. On the other hand, at very high temperatures, we expect 
a gas of almost free quarks and gluons. Results obtained in these limits, which are both 
far from the temperatures of deconfinement and chiral symmetry restoration, cannot help 
to explore the nature of the transition itself. 

For such an exploration much of the effort done so far is based on the mean field 
approximation to the effective potential at finite temperature for four-fermion models or 
other low-energy effective models of QCD. Unfortunately, the mean field approximation 
neglects fiuctuations, and it misses completely, in such a way, the role of bosonic contri- 
butions to the thermodynamical potential. These however should be, as already said, the 
relevant ones at very low temperatures. 

Recently there have been attempts to recover some information missing in the mean 
field approximation, by means of the expansion, where labels the number of fermion 
species 0, ^. In fact the role of scalar and pseudoscalar fiuctuations appears at the first 
order of the expansion and allows to test their relevance in various regions of temperature. 
For the moment being these attempts have been made for the Nambu-Jona Lasinio (NJL) 
model 0. 

Here we present an application, based on a general scheme which allows for the 
expansion of the effective action in four-fermion models [§], to the Gross- Neveu model in 
1+1 dimensions p. We include in the Lagrangian a current fermion mass to avoid infrared 
divergencies, and include the one pion-like composite states and the scalar composite. To 
verify whether their contribution to the free energy is dominant at low temperatures, we 
have to evaluate the effective potential at least at order 1/A^. In the imaginary time 
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formalism, we have to sum up over discrete Matsubara frequencies twice: the first time 
for fermions, the second one for bosons. 

The effective potential at this order of the expansion appears as a sum of various terms. 
The first is the zeroth order, mean field term. It has a zero temperature part and a finite 
temperature part, obtained by summing over fermionic energies. 

A second term is a pure temperature term, obtained by summing over both fermion 
and boson energies. It contains a Bose distribution function and it is expected to dominate 
in the expression for the pressure at low temperatures. In the following we show this to be 
indeed the case, and that in particular the largest contribution comes from the integration 
over almost on shell external pion momenta (close to the pion propagator pole). 

The last term is a mixed one, a zero temperature term with respect to the sum over 
Bose frequencies, depending on temperature through the bosonic propagator, obtained by 
summing over Fermi frequencies. 

This part of the effective potential presents non trivial aspects which are interesting 
as far as the model itself is concerned, but its detailed study would go beyond the scope 
of the present work. Actually we are primarily interested in showing how, by applying a 
general formula for four fermion theories one can evaluate the pressure at order 
and exhibit the role of pions at low temperatures. Since the mixed term is subleading for 
very low temperatures and small quark masses (see sect. 5), it does not appear in our final 
results. 

In the following paragraph we give a review of the expansion of the partition 
function, for the case of the Gross-Neveu model in 1+1 dimensions. In sect. 3 and sect. 4 
we extend the formalism to finite temperature, and calculate the pressure at order 1/A^. 
In sect. 5 we show the temperature behaviour of the pressure for the pure bosonic term 
discussed above and compare it to the mean field term. The results are summarized in 
the conclusions. 

Finally we summarize some useful calculations in the appendices, namely the ana- 
lytic properties of the inverse bosonic propagator (App.A), and the method employed for 
summing over the bosonic frequencies (App.B). 



2. expansion 



We start from the lagrangian of the massive Gross-Neveu model in 1 + 1 dimensions 




(2.1) 
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where is a multiplet of N degenerate fermion fields, and M is a bare fermion mass which 
exphcitly breaks the chiral invariance. 

The expansion can be carried out as follows 0. According to standard methods 
the generating functional can be rewritten by introducing a scalar field a and a pseu- 
doscalar field vr 0. At vanishing sources, we have 

1 



i d X 



^{id - M + ga + ig-i-,'n)i) - ^(a^ + tt^) 



(2.2) 



By carrying out the integral over the fermion fields, we obtain an effective action for 
scalar and pseudoscalar fields 



V{a)V{Tx) e'^B ^ 



Now we shift the cx-field, ga ^ ga + M, and write 



d X 



Ma 
9 



— i log det{id + ga + ig'-i^n) 



(2.3) 



(2.4) 



where 5Z is a counterterm which renormalizes the fermion loop [Q. 
Finally we redefine the parameters of the model as follows 

M 



X = Ng^ ; a = 



Ng' 



= {ga,gn); a = {l,i-f5) 



(2.5) 



Thus, the lagrangian for the 0-field is given by 



d X 



2A 



;i + 6Z) - Na<pi 



- iNTi \og{id + <p-a)= / d^xCB{,<P) (2.6) 



where we have used the identity log det = tr log, and carried out the trace over fermion 
indices. 

Up to now we have simply integrated over the fermion fields, passing from a description 
in terms of elementary fermionic constituents to a description in terms of bosonic com- 
posites. At this point we are ready to expand the generating functional Z in eq.( p.3| ) in 
series of 1/A^. For this we follow the method of ref.|T^, that we briefly review (notice that 
we are now performing an 1/N expansion and not the formal ^-expansion of the original 
reference). 

Let us couple an external constant source J to the field (p{x), and expand (j){x) around 
its vacuum expectation value (f 



(j){x) = (j)'{x) + if 

^[J] = (O+|0(x)|O")j 



5W 
Jj 



(2.7) 
(2.8) 
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Since J is constant, from the Lorentz invariance of the vacuum, it follows that is 
constant too. By supposing (^[J] to be invertible, we can define the effective action T[0\ 
as the Legendre transform of W[J] 

r[^] = W[J[^]] - j dPxJ[i{)]^ = -V(^) J dPx (2.9) 

From the definitions it follows that, at J = 0, must satisfy the stationary equation 

= (2.10) 



Sip 



At this point W[J] can be separated into two terms 

W[J] =Wo[J] + Wi[J] 

with iyo[<^] given by 

Wo[J] = J (fx[CB{^) + J^) 

and Wi [J] solution of the integral equation ||10| 



Wi[J] = -Hog J V{(j)')exp i J 



6(p 



(2.11) 
(2.12) 

(2.13) 



(2) 

The lagrangian {(f)',ip) is Cb{4>' + 'p) minus the constant and linear terms in 0', and 
can be obtained from eq. (|2.6|) . We have 

Sb = j d^x - — {(t)' + ^f{l + 5Z)-Na{(l)\ + ^i) - iNTilogiid + ^■a + cj)' -a) (2.14) 
Furthermore, by defining 

{id + ^- a)G{xi - X2) = 5^{xi - X2) (2.15) 

it follows 

Tr log{id + ip-a + (f)' -a) = Tt \og{id + ■ a) + Tr log {1 + G(f) ■ a) (2.16) 

Therefore, the standard mean field fermionic term in eg. ( p. 12 ), at J = 0, is 

N 



Wo 



d^x 
d^x 



~f{l + 5Z)-Na^i 
N 



iNTi \og{id + (p ■ a 
d^p 



2 2\ 

V ) 



(2.17) 



This is the only term which survives in the N +00 limit, since Wi is of order 1/N with 
respect to Wq. 
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The evaluation of Wi requires to expand Sb in eq.( p.l4| ) in powers of 0', starting from 
the quadratic terms. By rescahng the fields, (f)' (f)' / VN, and by ordering the terms in 
the exponent of eq.( |2.13|) in powers of it turns out that the leading term is just the 
quadratic one. 

Thus, by keeping only the leading term of the expansion, the field-dependent part 
of Wi can be obtained by 



exp <; -— / d XI 



1 + 6Z) 



— - j (fxi(fx2TT G{xi — X2)4>'{xi) ■ aG{x2 — Xi)(j)\x2) ■ a | (2.18) 

By going over to the Fourier space, we can diagonalize the exponent in the previous 
equation, thus obtaining a simple gaussian integral. 
We finally have 



J=a,iT 



(2.19) 



where iD^jlp) is the zero temperature inverse bosonic propagator, whose explicit form is 
given in Appendix A. 

At this point it is straightforward to verify that, apart from the linear term in <^i, 
both Wq and Wi depend on the chirally invariant combination = (7^(0"^ + vf^), which is 
the solution of eq.( |2.10| ). From this condition it follows that the pseudoscalar component 
vanishes, i.e. 

^={ga,0) (2.20) 



3. Finite temperature formalism 

The finite temperature extension is straightforward. For imaginary times we have 



J (fx — > J dx = -i(3V ; J 



E 



dp 
2^ 



(3.1) 



(277)2 ^ 

while the finite temperature effective potential, obtained as the generalization to finite 
temperatures of eq. (p.9|), is just the free energy density 

logZ 



f3V 



(3.2) 



According to the 1/iV expansion discussed in the previous paragraph, we can write 



(3.3) 



j=(T,7r 
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where is the purely fermionic standard mean field term coming from the one-loop 
calculation [ojn = {2n + 1)11/(3) 



= ^^^(1 + SZ) + Na^-^Y.J^ log + P' + -^l 



(3.4) 



and J-'^ is the contribution of bosonic fiuctuations {un = 2n'K / (3) 



B 



2/? 



E 



dp 
2^ 



log 



D- ^(iz/„,p) 



(3.5) 



The expression of the finite temperature inverse bosonic propagator iD^^ and its ana- 
lytic properties are given in Appendix A. 

Expression ( |3.5| ) and ( |3.4[ ) put in eg. ( p.3| ) give the equation of state, when evaluated 
at the solution of the gap equation ( p.lO| ) at finite temperature 

dH^^T) 



dip 







(3.6) 



which determines the evolution of the fermion condensate with temperature. 

By following a standard renormalization procedure and evaluating the sums in the 
last term of eq.(^T3D |TI[], the fermionic term can be written as 



|! (log 4 - l) + - A /*~,p,og ( 1 + 

47r \ mg / vrp h \ 



(3.7) 



where mo is the dynamical fermion mass generated at a = 0, T = 0. 
The sums that appear in the bosonic part of the free energy, eq. 



can be trans- 



formed into integrals by standard methods |T2[. In this procedure we must pay attention 
to the fact that the finite temperature inverse bosonic propagator iD~^{uj,p) has a cut 
across the origin in the complex cj-plane (for \uj\ < p), besides the cut extending from the 
continuum threshold to +oo (that is present even at T = 0) p, A detailed derivation 
is shown in Appendix B, from which it follows that the bosonic term can be cast in the 
form (uBiuj) = (e^^ - 1)"^ 



+ 00 



oo 277 Jo 
1 



+ 1™—/ IT dujnB{uj)\og 

e^O ZTII J-oo Ztt Jo 

apart from infinities independent of jS and cp. 



iDj ^{lu + ie,p) 
iDj^{uj — ie,p) 



(3. 
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Furthermore, since the function g{z) = iD ^{z^p) satisfies the Schwarz reflection prin- 
ciple, g*{z) = g{z*), by reintroducing the exphcit dependence on ip and T, we can write 

1 

^ dp 



+00 (~lp r+oo |- 

dulog iDj {iu!,p; ip,T) 



2ti Jo 



+ hm — 



2n Jo 



duj nsiuj) aig (iDj ^ {00 + ie,p;(p,T)j — TT (3.9) 



where arg{f) is the argument 6 G [0, 27r) of the complex number / = |/| exp(z6 



4. Free energy density and pressure at order 

According to the expansion discussed in the previous paragraph, can be written 
as the sum of the purely fermionic mean field term = NJ^o, and of the contribution of 
the bosonic fiuctuations ^ J-"^ = T\ 



j=(T,7r 



^((^(T),T) = iv^o(^(r),T) + j-i(^(r),T) 



(4.1) 



where V3(T) is the solution of the gap equation (|3.6| ). Therefore, in order to get a true 
expansion of the free energy density in eq. ([4.1|) , we need also to expand (^(T) in series of 

^(T)=^o(T) + ^^i(T) + --- (4.2) 

By requiring eq. (|3^) to be satisfied at each order in we obtain that <^o(7') is the 

solution of the mean field gap equation 

5-^0 (<^,T) 



dip 







(4.3) 



whereas </'i(T) is given by 



dip 



d'J'oi^,T) 



dip^ 



(4.4) 



Finally, by inserting eq.( [4.2| ) into eq.( [4.1| ), we have 

9^0 



^((^(T),T) = iV^o(^o(T),T) + 



dip 



(^i(T) + ^i((^o(T),T) + 

^o(r) 



= iV^o(<^o(T),T)+^i(<^o(T),T) + --- 
Let us now consider the pressure per fermion species, at order 



(4.5) 



V{T) 



-TmT),T)/N 



Po(^o(T),r) + -Pi(^o(r),r) + 



(4.6) 
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According to the previous paragraph the mean field pressure is given by 



^0(^0, T) 



,o2 



,o2 



A , log — - 1 - ay^o + 



pplogfl + e-^v/^) 



2 



(4.7) 



whereas the contribution at order can be separated into two terms 

1 



j=a,TT 



27T 



+00 dp r+00 
00 271 







(iculog iD - ^o,T) 



- ^ lim - / ^ I duj nsiuj) [arg (iL> • ^(t^ + ie,p; v^o, T)) - tt 



j=(T,7r 



TT J-oo 27r Jo 



In the next section we show the results for the temperature behaviour of the pressure, 
expanded as in eq.( |4.6| ), within the approximation of low temperatures and small current 
quark masses. 



5. Results 

Let us first comment on the second term in eq.( |4.8| ), which is the one expected to give 
the dominant bosonic behaviour at low temperatures. Actually such a term contains the 
Bose distribution function nsiuj), and, if iD^^ were the inverse free bosonic propagator, 
it would correspond to the standard pressure for a free boson gas 0, and as such it would 
contain the typical temperature behaviour of a gas of spinless particles. 

Even in the case of the inverse propagator having a more complex structure, the major 
contribution to the pressure at low temperatures should come from the integration over 
the region containing the poles corresponding to the composite bosons. In particular, 
since the pion corresponds to the lowest bosonic mass excitation because of its pseudo- 
goldstones nature, it characterizes the thermodynamics at low temperatures, where the 
leading behaviours are of the form ~ exp(— m/T), for the particles of mass m. 

Notice that in the present work we are including an explicit symmetry breaking, avoid- 
ing infrared divergencies related to the low dimensionality of the model. If one could 
go continuously to the massless case, one would recover the power law behavior in T of 
the thermodynamical quantities in the pion sector, evidencing the role of the pions as 
Goldstone bosons associated to the spontaneous breaking of the chiral symmetry. 

In conclusion the pion term in the second sum of eq. ( |4.8|) is expected to dominate the 
total pressure in eq. ([4.6| ) for low temperatures and small values of a = M/Ng'^, which is 
responsable for the explicit breaking of the chiral symmetry. In Figs.|l]-§ we show that this 
is indeed the case, in the present model. 
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We notice that the Gross-Neveu model is pecuhar since in presence of a bare fermion 
mass no bound states are present in the scalar sector (see Appendix A). However this 
particular feature of the model should not affect its physical interest to mimic some aspects 
of QCD at low temperatures, where the contribution of the scalar pole is depressed since 
m„ » m,r- 

Let us now consider the first term in the r.h.s. of eg. ( [4. 81) . Even this term has a 
temperature dependence, this time not through the Bose distribution function, but through 
the Fermi distribution functions appearing in the bosonic propagator at finite T, obtained 
by summing over fermions frequencies. 

This term can be separated in a zero temperature part 

^? = - E i / ^ r^du; log \zV,]{zu;,p;^o)] (5.1) 



27T J 27r 

J=CT,TV 







and a finite temperature part 

^ 2n J 271 Jo ^ 



j=o-,7r 



-r)-i 



(5.2) 



with an analogous separation for the inverse bosonic propagator (see eq.( |A.3D 



Vr^ ^ Do:] + -^^i (5-3) 

Vi is minus the zero temperature effective potential of order 1/A^ (evaluated at the 
minimum (fo). It is UV divergent and needs to be renormalized. The renormalization 
procedure for the massless Gross-Neveu model has been studied in literature [0, 0, P^ , 
although an explicitly invariant expression for the renormalized effective potential at this 
order has not been givenQ. 

However this term (once renormalized) is subleading, at low temperatures, with respect 
to pion pole. This can be verified by performing a low temperature expansion, and taking 
into account that depends on temperature only through (Pq{T), which has the following 
behaviour for low temperatures 

<^o(r) = (^o(o) + 5<^o(r) + ... (5.4) 

with 

^o(O) 



5^o{T) ~ VTe T (5.5) 



"'^We remind that since a bare fermion mass is present, there are no infrared divergencies in the pion 
sector, that otherwise would be present p6|. 
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Thus, as T goes to zero, Vi is exponentially depressed with respect to the pressure term 
coming from the integration around the pion pole, since <^o(0) > Tn^^. 

Let us now consider . We have verified numerically that at low temperature (and 
small enough values of a) also this term is subleading with respect to the pion pole term. 
This can be in part understood from the fact that its temperature dependence comes from 
the Fermi distribution functions, that, again, have an exponentially decreasing behavior 
at low temperatures, characterized by <^o(0)- 

In addition both and "Pf exhibit an imaginary part below a certain value of the 
quark condensate. This problem has been already noticed in literature for the zero temper- 
ature effective potential and we will not go further in its discussion. A complete 
calculation goes beyond the scope of this paper, since the problems are specific to the 
Gross-Neveu model, whereas the dominance of the pion pole is surely a general feature of 
QCD. 

We now discuss the behaviour of the pure temperature pion pressure term discussed 
above (see second line of eq.(^l8|)) 

\ r+QO r+co r / 1 \ 1 

'Pi {'fo,T) = -\im- — duj UBiuj) aigl iD:^ {uj + ie,p;ipo,T)) - IT (5.6) 

e^O vr J-oo In Jo L \ /J 

and we compare it to the mean field term, eq. ( |4.7| ). The results are shown in Figs.|I]-^ 
where we plot various terms of the pressure (divided by the square temperature) vs T/niQ 
for different values of the chiral symmetry breaking parameter a, and of the number of 
fermion species. 

According to the previous paragraph, and by subtracting the bag constant 7^o(<^o(0), 0), 
the expression for the pressure at order 1/A^ that we consider is 

ViT) ^ Vo{^o,T) - Po(^o(0),0) + ^Vn^o,T) (5.7) 

Furthermore, to better identify the role of the pion pole, we have divided the in- 
tegrations that appear in eq. ( |5.6| ) into two regions: (i) p < lu < Jp^ + A^pQ and (u) 



{uj"^ — p^){uj'^ — p^ — 4(^g) > 0. The region (z) is the one which contains the pion pole, and 
it is from here that the pressure takes its major contribution at low temperatures, always 
for small values of a. 

In Fig.|l](a) we show, for a = O.Olmo, the contribution to the pressure of the pion term 
Vi in eq. (|5.6| ), decomposed as follows: contribution of the pion pole (continuous) and that 
coming from the integration over the region (u) (dashed). In Fig.^(a) are shown the same 
quantities of Fig.|^(a) , for a = O.OOlmo. 

As anticipated in the previous paragraphs, we see that at low temperatures the pseu- 
doscalar bound state (provided the current quark mass is small enough) dominates over 
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the mean feald pressure, reflecting the fact that the pion is the lowest mass state by virtue 
of its nature of pseudogoldstone. Its contribution grows more and more as we decrease the 
bare fermion mass. In the massless case we would recover a power law behaviour in the 
limit T — s> 0, instead of the exponential one shown in Figs.|l]-0 (i.e. the curves ofV/T"^ and 
the one of VJ /T"^ containing the contribution of the pion pole would tend to a constant 
value and not to zero). 

In Fig.|I](6) we plot V/T"^ (see eq. (|5.7| )) for = 3 (short-dashed), = 10 (medium- 
dashed), and in the mean field case = +oo (continuum line). In this figure a = O.OlmQ. 
The analog results for a = O.OOlmo are shown in Figs.|^(6). 

Finally, in Fig.|^, we compare to the sigma contribution coming from the second 
term of eq.(|4.8|) for a = O.Olmo (a) and a = O.OOlmo (&). Although, as already noticed, 
the sigma term is peculiar in this model (since there is no scalar bound state) we show this 
figures to make manifest that its contribution is negligible at low temperatures, whereas, 
close above the critical temperature Tc = 0.57mo, the scalar and pseudoscalar fields have 
to be practically degenerate, since the chiral symmetry is almost restored (there is no real 
restoration due to the explicit symmetry breaking term). 

6. Conclusions 

We have applied a general formula, valid for four-fermion theories, to calculate the 
finite temperature effective potential at order 1/A^ in the massive Gross- Neveu model in 
D = 1 + 1. We have shown that the pion-like composite dominates the thermodynamics at 
low temperatures, correcting significantly the mean field approximation. We have shown 
in particular that the leading contribution at low temperatures comes from integration 
close to the pion pole. This feature, which is expected to be also general to QCD and 
related models, emerges clearly from our calculations. 



This work has been carried out within the programm Human Capital and Mobihty 
(BBW/OFES 95.0200; CHRXCT 94-0579) 
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Appendix 



A. Analytic structure of the bosonic propagator 

The zero temperature inverse propagator in eg. ( p.l9| ) is given by (we omit the index j) 



(A.l) 



[l + 5Z)-2i\ j 



(27r)2 [(g^ + p^) (g/^ + p/^) - {q^qt^ - cp^) \ 

where ± refer to the scalar and pseudoscalar self energy Il{uj,p) respectively, and 
{uj,p). By renormalizing in a standard way we have 

'^iDQ^{u,p; Lf) = -log^ 
A rriQ 



log I. . , . . — — :- (A.2) 



'(Cj2 - p2) (cj2 _p2_ 4^2) ^(^2 _ ^2) (^2 _ ^2 _ 4^2) + (^2 _ ^2) 

where e|,j^ = 4(^2 and e\j^ = 0. 

At finite temperature the inverse propagator D^^ appearing in eg. (|3r5|) can be calcu- 



lated by usual methods 



dq np{Eq 



27r E„ 



( I ■ \2 2 2 

-^2 V^^^"'' -V - Hi 



+ 



1 



1 



+ 



1 



(E, - zz/„)2 - [E, + iv^Y - El^^ 

1 



where 



{E^^^-iv^f -El {Eq+p + iPn? - El_ 

1 



Un = 2n7r/f3 ; upiE) 



E„ 



g2 -|- (^2 



(A.3) 
(A.4) 



Eg.( |A.3| ) can be cast in a more compact form 



-iD'^{iUn,p;'f,T) = —iDQ\iUn,p;(f)-4: / dq 

A JO 



E, 



(A.5) 



In order to evaluate the second term of eg. (|4.8| ) we take the inverse bosonic propagator 
as a function of u + ie instead of iun, namely iD~^{uj + ie^p; ^p,T). Then, the real and 
imaginary parts of iD~^ can be evaluated, in the limit e ^ 0, by using 



1 1 

lim = PV iirdix) 

e^o+ X + ie X 



(A.6) 
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Defining A = (o;^ — — A^p^) / (o;^ — p^), we have 



Re 



-\og^- arctan [^^^^ - + A l({uj, p) A < 



rrin 



log^ + VAlog 



x/A- 1 



-7/^ + A/f(a;,p) 



(A.7) 
A > 



Im 



27r ^ 1 



ro 



sinh f 

\ 2 , 



cosh ( ^ j + cosh ^ 



A < 



A > 



(A.8) 



27r ^ 1 



mn V-A 



arctan 



V=A 



/'^ + /f(w,p) A<0 



Va- 1 



VA + 1 



-/^ + /f(cu,p) A>0 



(A.9) 



Im 



ro 



— < 



TT 



sinh f 

I 2 , 



cosh I — 1 + cosh I — ^ — 



where we have defined the following integrals 



7^ = 4 / dq- 
Jo 



A < 



A > 



(A.IO) 



(A.ll) 
(A.12) 



tin 



(a;2_p2)2_4^2^2_4^2^2 



^9 ((a;2_p2)2_4^2^2 + 4p2^2^ - lepVl^^^-P^)^ 
When A > (and a;^ > 0) the last integral can be rewritten as follows 



7f = j= PV / (ig 

4a;p\/A -'o 



/o K 



(A.13) 



{{u^-p^f-Au'El-Ap-q^) (^-^ 



-2 9 -2 

T+ r - T- 
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with q±= {p± iuVAj /2 

Let us now summarize the main features of iD~^{z,p), with z = uj + ie, which follow 
fromeqs. (|X7D-(|CT|) 

a) at T = the inverse propagator iDq^ has a cut along the real z-axis, extending from 
the continuum threshold |a;| = y/p'^ + Aip'^ to +00 

b) for T 7^ a cut is present even for < p. This is due to the fact that at finite 
temperature "scattering" processes are possible between virtual scalar/pseudoscalar 
particles with space like momenta (cu^ — < 0) and quarks of the medium ^, [l^ 

c) in the pseudoscalar sector iD~^ has a zero, at least for low temperatures, for uj'^—p^ ~ 
amo, corresponding to the pion bound state [|T^ 

d) no scalar bound states are present for a 7^ 0, both at zero |19| and finite tempera- 
ture |1| 

Finally, we give the expression for the inverse propagator evaluated at imaginary en- 
ergies, which appears in the first term eq. ([4.8|) (A^; = (cu^ + p^ + Aip'^)l{ijj'^ + P^)) 



A 

2^ 



iD^'^{iuj,p]ip,T) = - log ^ + J A^; lo 



^-1 
/A^+1 



-I^ + AEl?ituj,plA.U) 



I^ + I?ituj,p) (A.15) 



B. Sum over bosonic frequencies 

The bosonic sums in eq.( |3.5| ) can be transformed into integrals by means of a standard 
procedure (for a review of finite temperature methods see JT^). Here we have to pay 
attention to the fact that, at finite temperature, the inverse bosonic propagator iD~^{z,p) 
(Re z = uj) has a cut extending along the real 2;-axis for —p < uj < p. 

Let us start by considering the following sum 

i^/(^z/0; yn = 2n'K/(3 (B.l) 

P n 

where we suppose that f{z) has no singularities on the complex plane apart from cuts along 
the real axis. Notice that in our case iD~^{z,p), although it has a cut for —p < u < p, is 
well defined in z = 0, since the discontinuity across the cut goes to zero as ut goes to zero. 



A 

2^ 



iD^^{iuj,p;ip,T) 



r 1 

log — + ^= log 



rrin 



'0 JAe 

for the scalars and pseudoscalars respectively. 



^+1 
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The first step is then to use the residue theorem to trade tlie sum for an integral over 
a circuit around the imaginary axis (see Fig.^. Since the integral along the whole contour 
r in Fig.^(a) gives not only the sum of the residues in z = iun, but also the discontinuity 
across the cut, we have to subtract it, namely subtract the contribution of the contour 
C = Ci + Cr (horizontal lines). 

Thus we can write 



P n 



1 

2TTi 



dzf{z)nB{z) - / dzf{z)nB{z 



c 



lim 



-ie+S 



+ioo+(5 



dzf{z)nB{z)+ dzf{z)nB{z) 

joo+<5 J +it+5 



+ lim 



dzf{z)nB{z) + 

+ioo—5 



where UBiz) is the Bose distribution function 



-too— (5 

dzf{z)nB{z) 



(B.2) 



nB[z) 



exp(/?z) 



1 



(B.3) 



which has poles in 2; = with residue equal to 

The third and fourth integrals of eq. ( [B.2|) can be transformed by a change of variable, 
z — > —z. Thus we obtain 



/5 



\ r-ie+S 

iVr.) = lim — / dz nB{z)[f{z) + f{-z)] 

5,6— >0 ZTTI J -ioo+5 



1 r+ioo+d 1 ^+ioo 

+ lim — / dznB{z)[f{z) + f{-z)] + — dzf{z) (B.4) 

where we have taken into account the fact that no singularities are encountered in per- 
forming the integral of f{z) on the imaginary axis (last term of eq.( p.4|) ). 

Finally, provided f{z)exp{—i3\z\) vanishes sufficiently fast at infinity and f{z) has no 
singularities on the complex plane (apart from the already mentioned cuts along the real 
axis), we can rotate the first integrals in eq. (|B.4|) along the circuit shown in Fig .^(6), and 
write 

1 1 r+oo '[ r+oo 

-73 fi^^ri) = lim — / du UbIuj) [f{uj + ie) - f{uj - ie)] + — / duf{iuj) (B.5) 
p „ f-*o m Jo ivr J-oo 



where we have also supposed f{z) = f{—z), and the limit 6 —>■ has been performed. 
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Figures captions 



Fig.l (a) Plot of the contribution to the pressure of the pion term in eg. ( ^.61) , decom- 
posed as follows: contribution of the pion pole (continuous) and that coming from 
the integration over the region (u;^ — p'^){uj'^ —p^ — 4<^q) > (dashed). These curves 
are obtained for a = O.OlmQ. 

(6) Plot of V/T"^ (see eg .15.71) for = 3 (short-dashed), = 10 (medium-dashed), 
and in the mean field case = +oo (continuous line). In this figure a = O.Olmo- 



Fig. 2 (a) Plot of the contribution to the pressure of the pion term Vi in eg.(5^), decom- 
posed as follows: contribution of the pion pole (continuous) and that coming from 
the integration over the region (u;^ — p'^){u'^ —p^ — 4<^q) > (dashed). These curves 
are obtained for a = O.OOlmo. 

(6) Plot of V/T"^ (see eg.|5.71) for = 3 (short-dashed), A^ = 10 (medium-dashed). 



and in the mean field case A^ = +oo (continuous line). In this figure a = O.OOlmg. 



Fig. 3 (a) Plot of the contribution to the pressure Vi in eg.(^]8|) of the pion (continuous) 
and of the sigma (dashed). These curves are obtained for a = O.OlmQ. 

(6) Plot of the contribution to the pressure Vi in eg. (^.8]) of the pion (continuous) 
and of the sigma (dashed). These curves are obtained for a = O.OOlmg. 

Fig. 4 (a) Contour integration used in eg.( [B.2| ). F labels the whole contour, while C = 
Ci + Cr is the horizontal contour circumscribing the cut across the origin. 

(6) Contour integration used in eg.( |B.4| ). F' labels the vertical contour along right 
side of imaginary axis, which is transformed in the horizontal contours Cr and 
circumscribing the cuts along the real positive axis. 
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